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TECHNICAL MEMORANDUM X-64605 


A FIRST-ORDER THEORY FOR ROTATING 
GRAVITY GRAD I OMETERS 

SUMMARY 


The equations of motion for a rotating gradiometer are developed and 
specialized to a forced harmonic oscillator model that may be solved analyti- 
cally. Numerical results are given comparing the analytical approximation 
to numerical integration of the exact equations. 


INTRODUCTION 


A rotating gradiometer may be visualized as two perpendicular 
dumbells coupled at their centers by a mechanical torsional spring pivot. 

G radiometers have been studied experimentally and theoretically for appli- 
cation to a variety of problems involving the measurement of gravitational 
fields [1-3]. 

Since the gravitational torque on each dumbell of the rotating gradi- 
ometer is opposite (and equal if they have equal moments of inertia) , there 
will be a deflection of the spring pivot coupling. This deflection is the 
quantity that is observed and related to the torque which may then be related 
to the external gravitational field and mass. 

The general relationship for the torque on any rigid body because of 
an external gravitational field is presented in the Appendix. This report 
begins with that basic relationship and develops the specific equations for a 
rotating gradiometer. Conditions at resonance are also discussed since the 
gradiometer may be rotated at a speed commensurate with the fundamental 
frequency of the spring to give a resonant condition. 



GENERAL EQUATIONS 


As shown in the Appendix, the gravitational torque on any rigid body 

is given by 


N =-G J V,(^Idm , (1) 

m 

where G is the universal gravitational constant, m is the external mass, 

I is the moment of inertia matrix of the rigid body, and V, is the gravity- 
eradient matrix of the gravitational field generated by m. It is also shown in the 
i t rgx that equation ( 1) is invariant under orthonormal coordinate 

rniations . 


if m is considered a spherical homogeneous mass and N, V, and 
are referred to the principal axes of the small body, equation ( 1) becomes 


N '= 


3 m G 

R 


X f R 


( 2 ) 


’ aj fie primes indicate that the quantities are referred to the principal 

as ol the small body and l' is the diagonalized moment of inertia matrix of 
’0 small body; i.e. , 


AGO 
0 B 0 

0 0 C 


Figure 1 gives the geometrical relationship of the prime system 
( fixed along the principal axes of dumbell no. 1), the unprimed system (fixed 
in the laboratory) , and the R vector. 

The torque on dumbell no. 1 referred to the x,y, z system may be 

written as 


N 


3mG 
fs- 


R'’ 


R X T I 


T — 
T R 


(3) 




Figure 1. Geometry of principal coordinates and laboratory coordinates. 

where T is the transformation matrix from the primed ( principal) axes of 
dumbell no. 1 ; i. e. , 

cos - sin 0 

T = sin Qi cos 0^ 0 

0 0 0 _ 

and 

'R cos d cos (p~\ R 

X 

R = R cos 9 sin (±> = R 

y 

R sin o' R 

^ z 



After denoting the principal moments of inertia of dumbell no. 1 as 
Aj, B|, andCj, it is assumed that A^= 0 and = Cj. Thus, ij becomes 


ro 0 0 

= p Cj 0 
0 0 Cj 


Using these results in equation (3) gives 


Ni = 


N 


( 1 ) 


X 


N 


N 


( 1 ) 


(i) 


3mGCi ^ . 0 ^ ^ ^ 

— n— ^ (R R sim0i + RR sin cos 
R y z ^ X z ^ '■ 


(-R R cos^0i+RR sin 01 cos 9 1 ) 
R X z ^ y z ^ ^ 


3m GC 


R' 


( R^ - R^ ) p sin 201 + RR cos 20- 
. y X 2 X y 


(4) 


Similarly, the gravity-gradient torque on dumbell no. 2 is 


N, 


to 

I 

■ 

I X 


N<^> 


j y 


to 


L ^ 

L 


— (R R sin^02 + R R sin 02 cos 02) 
Jbc y z X z 


Qyy% 

— (-R R cos^ 02 + R R sin 02 cos Oo) 
X z ^ y z ^ 


(5) 


SmGCa 


i(R^ - R^ ) — sin 2 02 + R R cos 2 02 
y X 2 X y _ 


Now assuming that the deflection of one of the dumbells with respect 
to the other is restricted to the x, y plane that coincides with the x',y' 
plane, then only the N components of the torque contribute to this 
deflection. This allows equating the third component of the torque to the rate 
of change of angular momentum along the z axis. Before doing this, how- 
ever, the restoring and damping torques will be formulated . 
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The restoring torque on dumbell no. 2 is 
= -ka 


{ 6 ) 


and the danaping torque on dumbell no. 2 is 

r(D) 




h O' 


where k and h are assumed constant. Likewise, the restoring torque on 
dumbell no. 1 is 


nP^ = k a 


and the damping torque on dumbell no. 1 is 

,(D) 


N 


h a 


(9) 


The sum of all the torques on dumbell no. 1 may be equated to , 

and a similar operation may be performed for dumbell no. 2. Thus, 


^ , • , SmGCi 

= h q; + ka + 


I (R^ - sin2 01 + R R cos2 0i 

L y X 2 X y 


and 
C262 - 


, . . . 3mGC? 

ha - ka + — - jr-g 
K 


- R^ ) X sin 2 09 + R R cos 2 09 
y X 2 X y 


( 10 ) 


( 11 ) 


It is observed from Figure 1 that 


7T 




( 12 ) 


then 


02 - 01 = a 
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Using these results in equations ( 10) and ( 11) , a differential equation for a 
and 6 1 ' may be written as 


a 


h 


k k 


3mG 


Cl* %)“ ‘I'cUU (K2-E^)(sm29; -sin20,) 


3mG „ „ , 

+ — R^Ry(cos 2 02 “ cos 2 0-^) 


(13) 


and 


h 


3mG 


'1= — a a- — « + ^ 


( - R^ ) sin 2 0 1 + R R cos 2 0 1 

y X 2 ^ X y ^ 


(14) 


7T 

Replacing 02 with 0 ^ + a + — from equation (12) gives the 

Ci 

following differential equations for 'a and "d\ as functions of a and 0^: 


■a- 


Ha - K a + 


3mG 


2R 


r (R 


R ) ( - sin 2 0 1 cos 2 a 


cos 2 di sin 2 a 


- sin 2 0 ;i) + 


3m G 
Rs 


R R ( - cos 2 01 cos 2 a 
X y ^ 


+ sin 2 01 sin 2 a 


- cos 2 01 ) 


(15) 


and 


;• h - k 

= c‘i “ 


f'(R^ - R^) X 2 01 + R R cos 2 0i 
y X 2 . X y 


(16) 


where 


TT d h , k k 

H = •7T + HT" and K = — + — 

Cl C2 Cl C2 


These are the exact differential equations for a spherical homogeneous 
perturbing mass generating a torque on the coupled dumbells where the 
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restoring torque is proportional to the deflection (a) of the spring coupling 
and the damping torque is proportional to the angular velocity ( a) of the 
spring coupling. 


FORCED HARMONIC OSCILLATOR EQUATIONS 


Equations (15) and (16) are not forced harmonic equations for the 
deflection, and certain assumptions are necessary before they can be reduced 
to forced harmonic equations. First, the small angle approximation must be 
made for a so that cos 2 o; = 1 and sin 2 a = 2a = 0 . The second 
assumption that is necessary is to assume that R^, R , and R are constant 

A third assumption is that = cot where co = constant angular rate. 

If the above assumptions are made, then 

GL+Hd'+Ko'= - (R^ - R^ ) sin 2 cot - R R cos 2 cot ( 17) 

R y X R X y 

and 


01 = cot . (18) 

Since R= (R ,R ,R ) is assumed constant, it may be assumed that 
_ X y z 

R lies along the x axis without any loss of generality since this is equivalent 
to assuming that the x axis is along R initially; then, 

• 3niG 

a 4-H(i + Kq! = sin 2 cot , (19) 

xl 

which is the forced harmonic oscillator equation. 



SOLUTION OF THE FORCED HARMONIC 
OSCILLATOR EQUATION 


The solution to equations of the form of equation (19) are well known. 
An excellent discussion of the solutions may be found in Reference 5. 

H 

The solution to equation (19) for — < K may be written as 


CL = e (A| cos qt + A 2 sin qt) + — sin(2ct;t - j3) 

P 


( 20 ) 


where 


b = 


H_ 

2 ’ 



n^ = K 


3mG 

R^ 


J 


[n“^ - {2co)^]^ + (4bo;) 


tan^ 


4ba; 

n^ - ( 2co) 2 ’ 


and A I and A 2 are initial integration constants. 


One quantity that is usually of interest is the time constant r, which is 

— folding time. This is 


sometimes referred to as the relaxation time or the 


the time required for the first term (transient term) to reach — of its 

1 2 ^ 

original value . This occurs when br = 1; i.e., t= — = — which is seen 

b H 

to depend on the damping of the oscillator. When t — «>, the transient term 
goes to zero. This is the transient time that is seen to take several time 
constants to get down to one -hundredth or one -thousandth of the initial value . 
After the transient time, the oscillator is considered to be in steady state with 

amolitude 


3mG 

R^ 

|[n^-(2o;)2]2+ (4bu;)2|i/2 


a 

P 


( 21 ) 



( 26 ) 


or for n » — = b 


Q 


g „ BI. 
2b 2 


and from equation ( 24) , 


2nr 


which gives 


2f = 


2Q 


( 27 ) 


(28) 


This indicates that a high f and high Q are incompatible. 


COMPARISONS AND RESULTS 


The results of the forced harmonic solution and numerical integration 
of equations (15) and (16) are shown in Table 1. The initial conditions 
assumed were as follows: 

K = (4'ir)^ radians/ sec^ , 

H = 0. 404 radian/sec , 

and 


CO = 2tt radians/ sec 

As can be seen, the steady-state results agree through the third 
significant figure. From these results, it may be concluded that the forced 
harmonic oscillator approximation is sufficient for design purposes. It should 
be noted that the forced harmonic oscillator theory assumes a linear spring 
constant and damping coefficient. 
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Another point of interest is the driving frequency a; required to give 
maximum steady-state amplitude. This occurs when p is a minimum; thus, 
differentiating p with respect to 2oj and equating the result to zero gives the 

conditions 

(2o;f= n^ - 2b^ 

or 


p = 2b^y n^ - b^ ( 22) 

for maximum amplitude. These results may be used for determining the best 
operating conditions of a gradiometer. 

Savet [3] has suggested a figure of merit f for various sensors defined 
as the sensitivity u divided by the time constant squared The sensitivity 
is defined as the steady-state displacement amplitude divided by the amplitude 
of the forcing term. This figure of merit in the present discussion turns out 

to be 



Using equation (22) in the above. 


2br‘ 






2t 




'V 


(24) 


If the quality amplification factor . Q is defined as the energy stored 
in the spring divided by the energy dissipated per cycle at resonance and 
steady state, then 


Q- 


2ir 


i-e 


-47t/ W n^-- 


? 


1-e 


2ir 

-8iri 


(25) 


9 



TABLE 1. POINTS IN STEADY STATE FOE 1 c/s 


Numerical Integration 

Forced Harmonic Oscillation Solution 

Time (sec) 

Delta (radians 
X 10~®) 

Time (sec) 

Delta (radians 



X 10 ) 

30.0125 

-0.59040669 

30.0000 

-0. 59015672 

30.2625 

0.59056961 

30.2500 

0.59034407 

32.0125 

-0.59143408 

32.0000 

-0.59142200 

32.2625 

0.59152283 

32.2500 

0.59154729 

33.0125 

-0. 59174430 

33.0000 

-0.59188750 1 

33.2625 

0.59180415 

33.2500 

0.59198996 j 

42.025 

-0. 59316773 

42.0000 

-0.59363485 1 

42. 275 

1 

0.59323604 

42.2500 

0.5936516 1 

[ 

Steady State | 

29.5125 

-0.59004683 

30.0000 

-0.59015672 j 


Figure 2 shows the waveform for 20 sec, which indicates that the 
transient term lasts for about 4t since t is 5 sec for this case. 
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APPENDIX 

THE GRAVITY-GRADIENT TORQUE AS A CONTRACTED 

alternating product 

THE GRAVITATIONAL TORQUE OF ONE RIGID BODY ON ANOTHER 


To establish explicit invariant relationships between the monaent of 
inertia of a small body, the gravity gradient because of the presence of larger 
bodies, and the torque on the small body, it is necessary to define the moment 
of inertia and gravity gradient in such a way that their transformation properties 
are known. 


The AAoment of Inertia 


From Figur^-A, the moment of inertias usually are defined as 




f (p^ - p? ) dm for i = j 
m 


(A-1) 


where p = (pj, p2,P3) , i = {1,2,3} , and j s {1,2,3} . To check the 

transformation properties of I ., an orthonormal Cartesian coordinate 

transformation M is defined such that p = Mp ' ; i.e. , 


and 


P, 


M. ,P , 

ip’ p' 


P. 

3 


M. ,P , 

n T 


Then, equation (A- 1 ) becomes for i^j 
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I.. = - f M. P , M. P , dm’ - - r M. p ,P , M. , dm 
IJ ip' p' jq' q' ^ ip’ P’ q' iq' 


-M. ,/ r p ip , dm\ M, = M. ,I , , M. 

ip'l^ p q I jq ip' p'q' ]q' 


(Ml') M 

iq’ jq' 1] 


(A~2) 


T 

Thus, I = MI'M for i # j . For i = j, the integrals in equation (A-1) 

are all of the form f p.p. dm , which will also transform according to 

•^11 

m 

equation (A -2) . 


The Gravity Gradient 


The gravity gradient is usually defined as the second partial derivatives 
of a potential function, say V(xi,X 2 ,X 3 ) . Then, 



9^V 


9x. 9x. 
J 1 


(A-3) 


To check the transformation properties of V, _ an orthonormal Cartesian 
coordinate transformation M is defined such that x = Mx' . Thus, 


X. 

1 


M. 

iP 


X 


X , = M X. 
P P 1 1 


= M. , X , , and x , = 

jq q q 


X. 

J 


M X. 

q J 3 



Now 


V, 


ij 



tr ' 

V_ p 

c t 8x. 

P 1 


d dv 


M 


= 


a^v 


q p 


ax. ax , pi 

] p 


M M 


V, , >M M = M v; , , M 
pq qj p 1 pi p q qj 


M (V, M ) 

pi ' p j 


T 

T’ T* 

(M ) v;m 


.. = (mv;m ).. 

ij 1] 


Thus, ¥, = MV;M . (A -4) 

The Torque as the CAP Product of V, and I 

From Figure 1, the element of force d(dF) on the element of mass 

drn because of the element of mass dm * is 


(N’dF) = - Gdmdm * =-Gdmdm=i= f (d) 

d 


However, since d-= r + p and it is assumed that p « r, then the i 
component of f(d) is [4] 


(A-5) 

th 


f. (d) = f. (r + p) = 1 (r) + 
1 1 1 


ai _ 

p. s f.(r) + V, .. p. 
9x, ] 1 13 3 


(A-6) 


r = (Xi,X 2 ,X 3 ) , p = (Pi,P2.P3) = 


c J _ i 

j- = and V = — 

1 ax. d 


Where 



Now, substituting equation (A-6) into equation (A-5) , 


d(di5 = -Gdmdnx* f(rO + V, 


The torque about 0 because of the force of dm* on dm is 


d( dN) = p X d ( dF) = - Gdmdm * p x f ( r) + V, p ] 
= - Gdmdm *p xT(r) - Gdradni* ^x V, ^ 


Integrating first over m gives 


dN = - Gdm* J p x f ( r) dm - Gdm>:= J p x V, p dm 


m 


m 


The first integral on the right is zero because of the choice of ( 
of mass for m;. The cross product in equation (A -7) is 


p X V, p = 


P 2 V. 3 .P. - P3V,3.P. 

P3V.,.P. - P,V,3.P. 

P,V,3.P. P3V.^.P, 


J 


Expanding and integrating over M gives 


dN = - Gdm'* 






Ij ^2j " 2j ^Ij 


(A -7) 

as the center 
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Gdm.=:= 


(V,I)23 - (V,I)32 
(V,I) 31 - (V,I) 13 
(V,I)i2 - (V,I)2i 


(A -8) 


As a notation convenience, it may be observed that the result in equation 
(A -8) may be obtained by placing the column or row vectors of V, and I 
in a determinant and using the conventional cross product rule except that the 
dot product is used when multiplying the column vectors. Using the Contracted 
Alternating Product (CAP) symbol © to indicate this operation between V, 
and I, one has 





63 


V2 • h 

- V3 ■% 

Vi 

V2 

V3 

- 

V 3 -Ti 

- Vi ■% 

Jl 

h 

h . 


,Vi -T2 

- V2 -Ti_ 


(A-9) 


which is the same as equation (A-8) since V. = (V, . V, . , V, ) , 

1 11 1 ^ lo 

I. = (I , I I ) and both V, and I are symmetric. Thus, the torque 

i li 2i 3i 

given in equation (A-8) may be written as the CAP product of V, and I as 


dN - - Gdm=;> V, @ I 


or 


N = - G 



(A-10) 


This is considered to be desirable since it has the effect of separating or 
isolating those effects that depend on the field which is in V, and those that 
depend on the physical characteristics of the body which is in I. It also gives 
their relationship to the torque. 
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The Torque Referred to a Coordinate System in Which 

1 is Diagonal 


It may be shown that applying an orthonormal transformation to the 
CAP product of two matrices (or tensors) is the same as the CAP product of 
the transformed matrices (or tensors) . It has already been seen that V, and 
I transform according to equations (A-2) and A -4) ; then, if M is the trans- 
formation from N' to N, 

dN = MdN' = -Gdm*MVj ©I = -Gdm* MV; 

= -Gdm* V,@ I 

or 

N= - G f V,0 Idm* . (A-11) 


From the preceding discussion it is clear that to refer the torque N 
to a particular coordinate system, both V, and I must also be referred to 
that system. For example, I may be diagonalized by referring I to the 
principal axes of M. In this system, the CAP product becomes 



V,23(I'22 “ I'ss) 
V;23a'33 - I'll) 
V:i2dii -I'22). 


(A -12) 


where the primes indicate that the quantities are referred to the principal axes 
of m . From equation (A-12) , the torque in the primed system depends on the 
three cross gradients; however, if the torque is transformed by a general 
rotation of coordinate axes, each component of the torque depends on all five 
of the independent gradient elements as indicated by equation (A-9) . 
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' li r-iue of a Spherically Symmetric Homogeneous Body on a 

Small Distant Body 

If m* is spherically symmetric, then p= (0,0,0) in Figure 1-A so 
that r = R - p* = R; then, 

?(?) - f(R - P) = f(R) = -|s 

9R3 R® 



Substituting the above into equation (A -10) and integrating gives 

N = R' X I'E’ , (A -13) 

R 

where 

(El, R2= R 3 ) 

This is the same result as used in Reference 6. 

SUMMARY 


An explicit relationship has been developed between the gravity gradient, 
moment of inertia, and torque of one large rigid body on another small body. 
From this theory, conclusions can be drawn concerning observable quantities 
and derived information from the observables. Also various applications, 
measLiremeEt concepts, and analysis techniques can be developed. 
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